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Abstract 

In the previous work, it was shown that, in supersymmetric (matrix) discretized quantum me- 
chanics, inclusion of an external field twisting the boundary condition of fermions enables us 
to discuss spontaneous breaking of supersymmetry (SUSY) in the path-integral formalism in 
a well-defined way. In the present work, we continue investigating the same systems from the 
points of view of localization and Nicolai mapping. The localization is studied by changing of 
integration variables in the path integral, which is applicable whether or not SUSY is explicitly 
broken. We examine in detail how the integrand of the partition function with respect to the 
integral over the auxiliary field behaves as the auxiliary field vanishes, which clarifies a mecha- 
nism of the localization. In SUSY matrix models, we obtain a matrix-model generalization of the 
localization formula. In terms of eigenvalues of matrix variables, we observe that eigenvalues' 
dynamics is governed by balance of attractive force from the localization and repulsive force 
from the Vandermonde determinant. The approach of the Nicolai mapping works even in the 
presence of the external field. It enables us to compute the partition function of SUSY matrix 
models for finite N (N is the rank of matrices) with arbitrary superpotential at least in the 
leading nontrivial order of an expansion with respect to the small external field. We confirm 
the restoration of SUSY in the large- N limit of a SUSY matrix model with a double-well scalar 
potential observed in the previous work. 



1 Introduction 



Spontaneous breaking of supersymmetry (SUSY) is one of the most interesting phenomena 
in quantum field theory. Since in general SUSY cannot be broken by radiative corrections 
at the perturbative level, its spontaneous breaking requires understanding of nonpertur- 
bative aspects of quantum field theory [Tj. In particular, recent developments in nonper- 
turbative aspects of string theory heavily rely on the presence of SUSY, which is however 
lost in the standard model. Thus, in order to deduce predictions to the real world from 
string theory, it is indispensable and definitely important to investigate a mechanism of 
spontaneous SUSY breaking in a nonperturbative framework of strings. Since one of the 
most promising approaches of nonperturbative formulations of string theory is provided 
by large- N matrix models [21 El H] (N is the rank of matrix variables), it will be desirable 
to understand how SUSY can be spontaneously broken in the large-iV limit of simple 
matrix models as a first step. For example, in IIB matrix model [3] it has been suggested 
that the rotational SO (10) symmetry is spontaneously broken in the large- N limit [5]. 
This tempts us to expect that SUSY is also broken in the large- iV limit of this model. 
Analysis of SUSY breaking in simple matrix models would help us find a mechanism which 
is responsible for possible spontaneous SUSY breaking in nonperturbative string theory. 

For this purpose, it is desirable to treat systems in which spontaneous SUSY breaking 
takes place in the path-integral formalism, because matrix models are usually defined 
by the path integrals, namely integrals over matrix variables. In particular, IIB matrix 
model [3] defined in zero dimension can be formulated only by the path-integral formalism. 
Motivated by this, in the previous work [BJ, we constructed the path- integral formalism for 
SUSY (matrix) quantum mechanics on discretized Euclidean time t G {1, 2, • • • , T}, which 
includes cases that SUSY is spontaneously broken 0. It is formulated in a well-defined way, 
by introducing an external field, which explicitly breaks the SUSY, to twist the boundary 
condition of fermions in the Euclidean time direction. In this setup, we compute an order 
parameter of SUSY breaking such as the expectation value of an auxiliary field in the 
presence of the external field. If it remains nonvanishing after turning off the external 
field, it shows that SUSY is spontaneously broken because it implies that the effect of 
the infinitesimal external field we have introduced at the beginning remains. Here, it 
should be noticed that, if we are interested in the large- N limit, we have to take it before 
turning off the external field, which is reminiscent of the thermodynamic limit of the 
Ising model taken before turning off the magnetic field in detecting the spontaneous Z2 
breaking. In the formalism proposed in [6], for the expectation value of an auxiliary 
field, the external field plays the role of a regulator by which it is computed in a well- 
defined manner. In particular, if we take the periodic boundary condition for fermions, 
the partition function is essentially the Witten index [7] which vanishes when the SUSY 
is spontaneously broken. However, since the external field explicitly breaks the SUSY 
by a small amount, the partition function with the external field becomes nonzero, and 
the expectation value normalized by the partition function is well-defined. Moreover, we 

1 Notice that SUSY can be spontaneously broken in systems denned in less than one-dimension as 
discussed in [5J. Namely, an analog of the Mermin-Wagner-Coleman theorem does not hold for SUSY. 
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have seen that the expectation value of the auxiliary field is also well-defined in the limit 
turning off the external field due to cancellation of its dependence between the numerator 
and the denominator. This shows how the expectation value of the auxiliary field can 
have nonzero value in the path-integral formalism. 

In view of this, it is quite important to calculate the partition function in the presence 
of the external field in the path integral for systems which spontaneously break SUSY. 
Especially it would be better to calculate it in matrix models at finite N in order to 
observe breaking/restoration of SUSY in the large- N limit. In this paper, we address this 
problem by utilizing two methods: localization [8] and Nicolai mapping |9J. As for the 
localization, we make change of integration variables in the path integral, which is always 
possible whether or not the SUSY is explicitly broken (the external field is on or off). It 
is investigated in detail how the integrand of the partition function with respect to the 
integral over the auxiliary field behaves as the auxiliary field approaches to zero. It plays a 
crucial role to understand the localization from the change of variables. To our knowledge, 
this kind of investigation has not been found in the literature. In the case of discretized 
SUSY quantum mechanics with Q-SUSY preserved, it implies that the path integral 
receives contributions only from the fixed points of Q-transformation and reproduces 
known results for the localization formula. In particular, for the T = 1 case corresponding 
to the zero-dimensional model, the fixed points of Q-transformation are nothing but the 
critical points of superpotential, i.e. zeros of the first derivative of superpotential. In the 
case of SUSY matrix models, analogous localization formula can be obtained. However, in 
terms of eigenvalues of matrix variables, an interesting phenomenon occurs. Localization 
attracts the eigenvalues to the critical points of superpotential, while the square of the 
Vandermonde determinant arising from the measure factor prevents the eigenvalues from 
collapsing. The dynamics of the eigenvalues is governed by balance of attractive force 
from the localization and repulsive force from the Vandermonde determinant. Without 
the external field, contribution to the partition function from each eigenvalue distributed 
around some critical point is derived for a general superpotential. In the case of a double- 
well scalar potential, it leads to the statement (4.17) in [6] in the large- N limit. When the 
external field is turned on, computation is still possible, but we find that a method by the 
Nicolai mapping is more effective. Interestingly, it works for SUSY matrix models even in 
the presence of the external field which explicitly breaks SUSY. It enables us to calculate 
the partition function at least in the leading nontrivial order of an expansion with respect 
to the small external field for finite N. We can take the large- N limit of our result before 
turning off the external field and detect whether SUSY is spontaneously broken or not 
in the large- N limit. As a byproduct of the analysis, we give a clear argument for the 
restoration of SUSY in a SUSY matrix model with a double-well scalar potential at large 
N, which was observed in [5J. 

This paper is organized as follows. In the next section, we consider change of variables 
in the path integral for discretized SUSY quantum mechanics leading to localization. It 
is pointed out how it works by investigating the behavior of the integrand of the partition 
function as the auxiliary field becomes small. In section |3j a similar method is applied 
to SUSY matrix models, and a matrix-model generalization of the localization formula 
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is derived. In section HJ we make an expansion of the partition function with respect 
to a small external field and derive a formula for a general superpotential in the leading 
nontrivial order of the expansion. It is valid for arbitrary N. By applying it to the case of 
a double-well scalar potential, we confirm the restoration of SUSY in the large- N limit of 
this model discussed in [6]. We summarize the result so far and discuss future directions in 
section Details of localization in discretized SUSY quantum mechanics with T > 2 are 
discussed in appendix^ Finally, some computational details are presented in appendix [Bl 

2 Change of variables and localization in discretized 
SUSY quantum mechanics 

As discussed in [6], in order to discuss spontaneous SUSY breaking in the path-integral 
formalism of (discretized) SUSY quantum mechanics or SUSY matrix models, we intro- 
duce an external field to twist the boundary condition of fermions in the Euclidean time 
direction and observe whether an order parameter of SUSY breaking remains nonzero 
after turning off the external field. This motivates us to calculate the partition function 
in the presence of the external field. In the following, we consider systems of SUSY quan- 
tum mechanics on discretized Euclidean time t G {1,2,- •• , T}. As shown below, it is 
possible to introduce such an external field even in zero dimension (T — 1). Therefore, 
by considering the simplest zero- dimensional models, it is expected that we can extract 
some essential properties of the partition function in the presence of the external field 
without touching technical complexity. In this section, we consider change of variables 
which leads to the localization of contribution to the path integral and will be useful in 
the computation of the partition function. 

2.1 Introduction of external field 

We begin with SUSY quantum mechanics whose action is given by 



where the Euclidean time direction is compactified by (3, (') means the time derivative, 
and W'{4>) and W"((p) are the first and second derivatives of the superpotential W((f>) 
with respect to (p. In this paper, we focus on the case that W(<j)) is a polynomial of <fi. S 
is invariant under one-dimensional M = 2 SUSY transformations generated by Q and Q, 
which act on fields as 



S 




) 



(2.1) 



Q(j) = ip, Qip = 0, Qip = -iB, QB = 



(2.2) 



and 



Qcj) = - 1 p j Qijj = o, Qtp = -iB + 20, QB = 2iip. 



(2.3) 



3 



They satisfy the algebra 



Q 2 = Q 2 = 0, {Q,Q} = 2d t . (2.4) 
The invariance of S follows from its Q- or Q(Q-exactness: 

S = Q J dt${lB-(i + W*{4>j)} = QQf dt(^ + W^Y (2.5) 

The partition function is defined by 

Z = J VBV(j)V^V^e- s (2.6) 

with the path-integral measure normalized as 

J Vte-tf***®' = J VBe-ti*$ B ®' = 1, J Thl>V$ e~ $ dt ^ {t) = 1. (2.7) 

It is pointed out in [2] that if we take the periodic boundary condition for all fields, (12. 6p 
is equivalent to the Witten index [7], which vanishes when SUSY is spontaneously broken. 
It means that the expectation value normalized by the partition function is generally ill- 
defined in such a case. Since the vanishing partition function originates from cancellation 
between bosonic and fermionic states, we will introduce an external field which explicitly 
breaks the SUSY, in order to resolve the degeneracy and to fix a single vacuum in which 
the SUSY is broken. It is analogous to the magnetic field introduced in the Ising model 
in detecting the Z2 symmetry breaking. Let us modify the periodic boundary condition 
of the fermions to a twisted one as 

if>(t + /3) = e ia i[>(t), 4>(t + p) = e~ ia i;(t). (2.8) 

Here, the twist a corresponds to the external field. It is shown that in the presence of a, 
( 12. 6 p does not vanish and the normalized expectation value of the auxiliary field {B) is 
well-defined. Moreover, it turns out that {B) a does not depend on a and therefore a — > 
limit is also well-defined [6] . In this sense, a plays the role of a regulator by which we can 
calculate the expectation value of an order parameter of SUSY breaking unambiguously. 
Thus, the external field a provides a framework for discussing spontaneous SUSY breaking 
in the path-integral formalism. 

As a discretized version of ( 12. 5p and (12. ip . we consider 

T 

s = g^#)|^(t)-(0(t + i)-0(t) + w'(0(t)))} 



T r 



E 



t=i 



l -B{tf + iB(t) {<P(t + 1) - 0(t) + W'(<P(t))} 



+ m {^{t + 1) - m + w H w))W)} , (2-9) 
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which preserves Q-SUSY but breaks Q-SUSY by the discretization d. Let us express by 
S a the action (12. 9p under the twisted boundary condition 

0(T + 1) = 0(1), ^{T + 1) = e te V(l). (2.10) 

Namely, 0(T + 1) and if)(T + 1) appearing in ( 12. 9 p are understood to be replaced with 
0(1) and e ia il}(l), respectively. Then the partition function is defined as 

Z a = J f[ (dB(t) d(f>(t) #(t) d$(t)) e~ s <*. (2.11) 

We will fix the sign convention of integrals over Grassmann numbers as 

#(t) iP{t') = [ #(t) = V- (2-12) 



In the simplest case T — 1, the action and the partition function are expressed as 
S Q = l -B 2 + i W(0) + ij {e ia - 1 + W"{<P)) if>, 

Z a = — - / dBd(pdipdi!e- Sa . (2.13) 
27r y 

We see that the effect of the external field remains even in the zero- dimensional model 
and breaks the SUSY. 



2.2 Localization in T = 1 discretized SUSY quantum mechanics 

As a simple example of localization, we first discuss the T = 1 system (I2.13P under the 
periodic boundary condition (a = 0): 

S = l -B 2 + iBW'((f>) + $W"(<l>)il>, 

Z = — - / dBd<pdipdipe- So . (2.14) 

27T J 

Sq preserves the Af = 2 SUSY 

Q0 = ^, QV = 0, Qtp = -iB, QB = 0, (2.15) 

and 

Q0 = -^, Q^ = 0, Qip = ~iB, QB = 0, (2.16) 

which are reduction of (12.21) and (12. 3p to zero dimension. 

2 When T = 1 and all the variables obey the periodic boundary condition, the action is nothing but 
the dimensional reduction of (I2.1j) and invariant under both Q and Q as seen in the next subsection. 
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Let us consider the following field redefinition H (B, 0, ip, ip) — > (B, 0, ip, e): 

= + e^, $=-ieB. (2.17) 

Note that, from the SUSY transformation (12.151) . these can be rewritten as@ 

= + eQ0, ^ = + eg^- (2.18) 

It implies that and ip are expressed as the SUSY transformation from and ip = 
respectively, and that the SUSY transformation parameter e is regarded as a fermionic 
variable instead of xp. Then, from the SUSY invariance of So, 

S (B, 0, xp, $) = S (B + eQB, + eQ0, V + eQtfj, + eQ</0 

= S (B, 0, ip, $ = 0) = ^i? 2 + i W(0), (2.19) 

which is independent of e. (This expression can be directly derived by using = (p—iipip/B 
obtained from f)2.17p .) Furthermore, since the Jacobian associated with (I2.17P is computed 

as 

dB d4 d%l)d^p = 4 dB dj> dip de, (2.20) 
B 

and the B- integral in the partition function looks singular at B = 0, we can say that the 
change of variables (I2.17P is always possible for B ^ 0. However, notice that, if other 
-B-dependence than the Jacobian (12.201) arises which makes the 5-integral nonsingular at 
the origin, (I2.17P is possible even at B = 0. We will see such an example explicitly below. 

By using (I2.19p and (12.201) . we find that the path integral of the partition function Z 
given in (I2.14p is localized at B = 0. Namely, if we divide the integration region of B in 
( I2.14p into the vicinity of B = 0: {B | \B\ < e} and its complement {B | \B\ > e} as 

Z = Zf ] + Z , (2.21) 

where 

Z^=f dBE (B), Z = [ dBE (B), 

J\B\<e J\B\>e 

Ho(-B) = — - / c/0#^e" 5 ° (2.22) 

3 The argument leading to localization from a field redefinition is based on Chapter 9.3 in [10 . 
However, the auxiliary field is not introduced there and the treatment of path-integral measure seems 
somewhat incomplete. For instance, the second term of (9.35) in [10] does not vanish in general, contrary 
to the claim in [TP] , 

4 We can also consider another field redefinition which is expressed as Q transformation, and the 
argument proceeds similarly. 
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with < £ <C 1, then Z is shown to vanish due to the trivial e- integral after the above 
change of variables. 

On the other hand, for the purpose of examining whether (I2.17P works even in com- 
putation of Zq , we have to take account of other Independence than the one in (I2.20p 
and to observe if 5-integral still diverges or not. In order to see the behavior of Ho(-B) in 
the vicinity of B = 0, it is instructive to try the change of variables for \B\ < e. We have 

E (B) = -—l-e-^ B2 [°° d$e- iBW '® [ # de. (2.23) 
2tt B ,/_ 00 J 

When W'(<p) is a polynomial of degree p (p > 2): 

W'{4>) = 9 P r + Qp-i^' 1 + ---+90, (2.24) 

we rescale as 



(j)=\B\-p<p' (2.25) 
to extract \B [-dependence from the 0- integral : 

d e -iBW'(4>) = _L_ / d(j) > e -i SS n(B) gp<t> ' P ^ + 0(£ 1/ P) ] (2 26) 

) \B\ P J —DO 

Note that, since < £, the first term of Vy(</>) in (I2.24p becomes the most important 
after the rescaling (I2.25p . Hence, we see that S (-B) is singular as near the origin, 

and can be expressed as 

x [1 + O^ 1 ^)] . (2.28) 

The integrals of ?/> and e vanish, while the 5-integral is divergent. Since the expression 
(I2.28P is of indefinite form oo x 0, it is found that the change of variables (12.1 7p is not 
appropriate to compute Zq°\ (For p = 1 case, because the (^-integral in (I2.23P gives 5(B), 
it is clear that the B- integral in Z^ is divergent.) The indefinite form of Z^ under the 
change of variables (I2.17P can be understood to reflect that Z^ possibly takes a nonzero 
value if it is evaluated in a well-defined manner. 



The (//-integral is computed as 



/I (p : odd) 

_ 00 # ' Sin( ^ ,P)= |^2sin(^)r(l + i) ( P :even). ^ 
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Unnormalized expectation values For the unnormalized expectation values of B n 
(n > 1): 

(B n )' = — — [ dB dcfidipdip B n e~ s °, (2.29) 

27T J 

we use the same change of variables to have 

(B n )' = — [ dB d~4>di) dlB n - 1 e-^- lBW '^ . (2.30) 
2ir J 

In contrast to the case of the partition function, we will see that the change of variables 
( I2.17P is possible for any value of B in evaluating (B n ) . 

Notice that the B- integral is not singular at B = for n > 1. In fact, in the region 
\B\ < e, the i?-integral after the rescaling (I2.25P gives a finite value: 

/ dB B n_1 ~p e"^ 2 = —j (1 + 0(e 2 )) (2.31) 
Jo n- - 

for p > 2 case of W'(<p) in (12.241) . In p = 1 case, the 5-integral is clearly finite as 

J dB B n ~ l 5(B) = 5 nA . (2.32) 

It indicates that the change of variables ( I2.17P is allowed for any value of B. 
Thus, for all p > 1 in W'(4>), we obtain 

(B n )' = (n > 1) (2.33) 

from the trivial e- integral: J de = 0. 

( I2.33P implies that the localization to B = is realized in such a way that Ho(-B) is 
proportional to 5(5) (without derivatives of <5(-B)). It can be directly derived as follows. 
Applying the Nicolai mapping X = W'(<f)) to 

E (B) = — [ dcj)e-^ B2 - iBW 'M W"(<j)), (2.34) 

27T J 



we have 



1 /"°° 

H ( J B) = e-^ 2 — ft / dXe~ wx = $5(B), (2.35) 



where Jj is the mapping degree of X = W' ((/)). For W'(4>) given as a polynomial (I2.24p 
with the degree p, 

sgn(^p) for p: odd 
for p: even. 
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Localization to W{4>) = Since 

00 1 / 1 \ n 



e 2 



1 / i \ n 

n=0 ' ^ ' 



which follows from ( I2.33P for an arbitrary parameter u, the partition function can be 
expressed as 

Z = — I ' dBd<P d^J dtp e -^-iBW(4>)-jw»m_ (2>38) 



2vr „ 

Note that Z does not depend on the value of u. Let us take u > to perform the 
^-integration first. Then, 



V 27TU 



1 e-^^'W 2 (2.39) 



a/27TM 

In the limit u — > 0, the factor ^i— e'^^'^ 2 becomes <5(W // (0)), which directly leads to 
localization to the critical points of the superpotential satisfying W(<f>) = 0. 

In the case that the superpotential is a polynomial and its critical points are nonde- 
generate (i.e. W = 0, W" 7^ at the critical points), the limit u — > yields a well-known 
formula of the localization: 

/°° W"((h\ 
# 5{W \<j>)) W" '(0)= £ 1^^, (2.40) 

where the sum is taken over the critical points. Comparing (I2.40p with the 5-integral of 
f)2.35p . we obtain 

<t>:W'(<t>)=0 1 VYVI 

The same result can be obtained by the one-loop computation around the critical points. 
Let 0c be a critical point of W(<f>) and ip be a fluctuation around C : 

<l> = (f> c + y/u(p. (2.42) 

Then, (I2.39|) becomes 

Z„ = 7= / ^e-^"^)V^"(0 c ) + O( v ^), (2.43) 

where contribution around each of critical points has to be summed if W'{(f>) has two or 
more critical points. It is easy to see that the one-loop computation of ip reproduces the 
RHS of (12.401) . Note that higher loop contributions are 0(y/u) and negligible in the u — > 
limit. Thus, the one- loop computation around the critical points of the superpotential 
gives the exact answer of the partition function in all order of perturbation theory. 
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2.3 Localization in the presence of external field 

Next, we consider the system (I2.13P with the twisted boundary condition (a ^ 0). The 
field redefinition ( I2.17p changes S a to 

S a = X -B 2 + iBW'{$) - i(e ia - l)eBip, (2.44) 

which has e-dependence due to the twist a. We again separate the integration region of 
B as 



z a = + z a , 

Z<® = f dB~ a (B), Z a = I dB~ a (B), 

J-e J\B\>£ 



>\B\>e 

E a (B) = -— / d<pdipd^e- Sa . (2.45) 



2tt . 

For B^O, the change of variables (I2.17P leads E a (B) to[§ 



2tt B 



1 



= ( e - - 1) — / d4>e^ B WW ' W - (2.47) 

27T J 

E a {B) and thus Z Q do not vanish in general by the effect of the twist e ia — 1. This 
suggests that the localization is slightly violated by the twist. In an exceptional case of 
W'(4>) being linear, Z a vanishes, because the 0- integral gives S(B) whose support is out 
of the integration region of B. Thus, the localization persists in the presence of the twist 
a when W'((j)) is linear. 

On the other hand, the twisted partition function Z a is computed without using the 
change of variables ( I2.17P as 

Z a = ^ J dBd4>e^ B2 - iBW '^ (e ia -l + W"{(j>)) 



i e 2 



6 Note that the second line of (|2.47[) is not valid for B ~ 0. In the expansion of the last factor 
e i(e ia -i)eB^ = 1 + ^ e ia _ i^ B ^ in the first line of ^f]f , we should not drop the first term "1" although 

it yields vanishing Grassmann integrals. The reason is that the B-integral in Z^ is singular and that 
the total contribution to Za^ is of an indefinite form oo x which cannot be simply regarded as zero. 
Since the corresponding term is nothing but Z^ , we find 



4°) = Z<°> + (e ia -I) dB d^e-i B2 - ww '^. (2.46) 
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The second term represents the effect of the twist, which is the sum of Z a and the effect 
of the twist on Za^ (i.e. Z^ — Z^). Note that Z = Z^ from the localization seen 
in the previous subsection. The second term in (I2.48P of course tends to zero in the 
a limit, but notice that it becomes relevant when the SUSY is spontaneously broken, 
namely Zq = 0. Let us take a closer look at the effect of the twist on Z^ 1 which is the 
contribution from the vicinity of B = 0. When W'(<f>) is linear (W'(<p) = g±(j) + go), 

-i re poo 

Z W-Z® = (e ia -l)— / dBe- 1 * 32 d<Pe- w ^ + ^ 

2tt J_ £ J_ oq 

(2.49) 



0i 



is not zero even in the e — > limit, while a similar calculation tells us that Z a = 0. In 
contrast, when W'(4>) is a polynomial (I2.24p with the degree p > 2, we can show that 
Za^ — Z^ vanishes as e — > 0. Similarly to the argument in the previous subsection, 
rescaling (f> — > \B\~p <fi' yields 

Z (0) _ Z (0) = _ 1} J_ f dB J_ -\& f d(f) , e -is g n(B) 9p ^ N + Q^/rt 

27r J-e \B\p J 

= (e ia - 1) - ( r dB e'^ B2 ) !°° d<f>' cos(<?„0 ,p ) 

K \Jo Bp ) J-oo 



X 



[1 + 0{eV*)] ■ (2.51) 
Here the 5-integral is not singular, and we find 

z (o) _ z (o) = q ( £ i-j\ ^ ( e o) (2.52) 

for p > 2. Thus, we conclude that the effect of the external field in the partition function 
is irrelevant in the vicinity of B = except the case that W(4>) is linear. Note that, when 
W'((f)) is linear, the fermion determinant does not contain field variables, and that ~E a {B) 
is proportional to H (£>) as B. a (B) = e ~* +gl S (■£?). This explains the persistence of the 
localization under the twist. 



Unnormalized expectation values For the unnormalized expectation values of B n 
(n > 1): 

(B n )' a = I dBd(j)di)di)B n e- Sa 1 (2.53) 

ZTT J 



7 The (//-integral is calculated as 

^osfn^'H = ^ 9.rns I 

2^/ V P 



dcf>' cos(g p( f>' p ) = —V 2 cos ( ^- ) T ( 1 + - ) . (2.50) 
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the change of variables (12.1 7p leads to 



(B n )' Q 



1) — I dBd4>B n e-^ B2 - iBW '^. 
2tt 



(2.54) 



Note that, since the B-integral is not singular, the change of variables is always possible. 
Rewriting as 



(e ia - r 



1 



2,7T 



3$ 



e 2 



$2 



(2.55) 



we can express it by the integral of the Hermitian polynomials: 



d4>H n {W\ct>))e 



<s>=w(<t>) 



-WW 2 



[2.56) 



where the Hermitian polynomials are defined by 

HJx) = 



i 9 d n 
dx 1 



■e 2 . 



The expectation value of B n normalized by Z a : 

(B n ) a = (B% 



(2.57) 



(2.58) 



trivially vanishes as turning off a in the case that the degree p of W'(<p) (I2.24p is odd, in 
which the SUSY is not spontaneously broken (Z 7^ 0). However, for even p where the 
SUSY is broken (Z = 0), taking the ratio of (gag) and (I2"3B1) we obtain 



(B n ) a 



-i) n r_ oo oo #F n (W(0))e-5 



00 dcbe-^'M 2 

00 ~ 



(2.59) 



which can take a nontrivial value. Note that, because the factors (e ia — 1) appearing in the 
numerator and the denominator cancel each other, the value of (B n ) a is not dependent 
on a |. 

The argument so far presented for T = 1 can be extended to the case of general T. 
We put the discussion in appendix [XI 



3 Change of variables and localization in SUSY ma- 
trix models 

In this section, we discuss localization in SUSY matrix models, which yields some new 
features not seen in the previous section. Let us begin with a matrix-model analog of 
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(TO 



5 = Q£;jVtr#)Us(O-(# + l)-0(t) + ^W))) 



^JVtr 



-5(t) 2 + iB(t) + 1) - 0(t) + W(0(t))} 



(3.1) 



where all variables are N x N Hermitian matrices H. Under the periodic boundary condi- 
tion, this action is manifestly invariant under Q transformation defined in (12. 2p . 

We will focus on the simplest case T = 1 below. Under the twisted boundary condition 
(j2.10p . the action is 



Ntr 



-B 2 + iBW\<f>) + i) {e ia - 1) V> + il)QW\<i>) 



and the partition function is defined by 



(-1) 



N 2 



d B d 



d N ^d N ^) e-* 



(3.2) 



(3.3) 



where we fix the normalization of the measure as 



/5e" Mr ^ 2 » = 1, (-1) 

Explicitly, when W'((f)) is given as in (j2.24jl . ( 13 .2p becomes 

1 _ P „ , 

-B 2 + iBW'ty) + i> [e ia - 1) V + Yl 9k Y, $ ^ 



d N2 ^d N2 ^p) e - Ntt ^ = l. 



S n = Nti 



k-l 



uk-e-i 



k=l 1=0 



(3.4) 



(3.5) 



Notice the ordering of the matrices in the last term. We see that the effect of the external 
field again remains even after the reduction to zero dimension. When a = 0, S a= o 
is invariant under Q and Q given in ( I2.15P and (j2.16j) . both of which become broken 
explicitly in S a by introducing the external field a. 

Now let us discuss localization of the integration in Z a . Some aspects are analogous 
to the discretized SUSY quantum mechanics with T > 2 under the identification iV 2 = T 
from the viewpoint of systems possessing multi-degrees of freedom, while there are also 
interesting new phenomena specific to matrix models. We make a change of variables 



4> + eip, ip = ip — ieB, 



(3.6) 



8 From the viewpoint of spontaneous SUSY breaking in discretized noncritical superstrings, SUSY 
matrix quantum mechanics with a cubic superpotential and an Af = 1/2 version of the T = 1 case of 
(13. lj) are discussed in [IT] and [12] . respectively. 
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where in the second equation, ip satisfies 

Nti(B^) = 0, (3.7) 

namely, ip is orthogonal to B with respect to the inner product (Ai,A 2 ) = NtT(A\A 2 ). 
Let us take a basis of N x N Hermitian matrices {t a } (a = 1, ■ ■ • , iV 2 ) to be orthonormal 
with respect to the inner product: Ntr(t a t b ) = 5 a b- More explicitly, we take 

with Mb = \ \B\ \ = y/Nti(W) the norm of the matrix B. Notice that in the present case 
(N is general) ip is an N x N matrix and that e does not have enough degrees of freedom 
to parametrize the whole space of ?/>, which is in contrast with the N — 1 case ( 12. 17ft but 
analogous to the discretized quantum mechanics with T > 2 (IA.2j) . In fact, e is used to 
parametrize a single component of ip parallel to B. 
If we write (I3.3P as 

d N2 BE a (B), E a (B) = (-if 2 J d N2 <p (d N2 4jd N2 ^j e~ s «, (3.9) 

and consider the change of the variables in S a (5), B may be regarded as an external 

variable. The measure can be expressed by the measures associated with ip and e 

as 

d N2 iP = -? r ded N2 - 1 $, (3.10) 
Mb 

where d N2 ~ l ip is explicitly given by introducing the constraint (13. 7p as a delta-function: 
d N2 ~ l i = (-lf-'d^^f-^NtriB^)) 



M B v r V 

N 2 \ „ N 2 



if 2 - 1 (n#) TrT, Ba ^- ( 3 - n ^ 



.. Mb , 

va=l / a=l 

■?/5 a and B a are coefficients in the expansion of ip and -B by the basis {t a }: 

N 2 N 2 
a=l a=l 

( 03.1 Op and (13. lip are analogous to (lA.lOp and (1A.9|) in the discretized quantum mechanics 
with T > 2, respectively.) Notice that the measure on the RHS of (I3.10p depends on B. 
When B ^ 0, we can safely change the variables as in ( 13 .6p and in terms of them the 
action becomes 



S a = Ntr 
with Q<p = if). 



l -B 2 + iBW'($) + i {{e ia - 1)V + QW'{$f) - {e ia - l)ze£^ 



(3.13) 
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3.1 a = case 

Let us first consider the case of the periodic boundary condition (a = 0). Similarly to 
flA.5j) . S a =o does not depend on e as a consequence of its SUSY invariance, because (13.61) 
reads 

= + eQ0, $ = $ + eQ$. (3.14) 
Therefore, the contribution to the partition function from B ^ 

Z a=0 = [ d N2 BE a=0 (B) (0<£<1) (3.15) 

J\\B\\>e 

vanishes due to the integration over e according to (I3.10p . Namely, when a = 0, the path 
integral of the partition function (I3.3P is localized to B — 0. 

For the contribution to the partition function from the vicinity of B = 

Z { °1 = [ d N2 BE a=0 (B), (3.16) 

J\\B\\<e 

we can repeat the same argument as in the previous section. For instance, when W'(<j)) 
is given by (12.241) of degree p > 2, rescaling as 

1 = ^V, $ = (3.17) 



we obtain 

N 2 



7 (0) 



1 



2tt/ Wn j^+p 



dAfB^re-^ ] / cM B [ d N2 <f>' e"*^^' 



x />> [ded N2 -^' e-^^E&^'W^- 1 ] [i + 0( e i/p)] f (3.l8) 



where the measure of the 5-integral was expressed in terms of polar coordinates in R N 

as 

N2 JR* /i\t 
^ 2i? = n^=(^) Mf-'dMsdQB, (3.19) 

a=l * ^ ' 

and = represents a unit vector in R w2 . Since the e- integral vanishes while the 

integration of Mb becomes singular at the origin, Z^l takes an indefinite form (oo x 0). 
When W'(4>) is linear (p = 1), the ^-integrals in (13. 16p yield 

/ i\N* f (N* \ AT 2 

Vl^il/ yiiBii< e J Mb 



x / d N2 iP / ded^-^e"^^, (3.20) 



which is also of indefinite form - the I?-integrals diverge while j de trivially vanishes. 
Thus the change of variables ( 13 .6|) is not suitable to evaluate Z^Iq which possibly takes a 
nonzero value. 
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3.1.1 Unnormalized expectation values 

Next, let us consider the unnormalized expectation values of jjtrB n (n> 1): 

(^tr^y = J d N2 B (yriT) E a=0 (B). (3.21) 

Since contribution from the region | \B\\ > e is shown to be zero by the change of variables 
( 13. 6p . we focus on the B- integration around the origin (\\B\\ < e). 

When W'((j)) is a polynomial (I2.24p of degree p > 2, after the rescaling (I3.17P we 
obtain 

^tr^y = in*dtf B ArZ~ l ~'e-$»z\ Y N [l + Oie 1 

Y N = (^=) N Jdn B ^tr{Ql) J d^' e -*"M<W*) 

x f d N2 ij [ ded^-y e"^*^ 1 * - ^" 1 ]. (3.22) 



The A/"b- integral is finite, and it is shown that definitely vanishes in appendix [B] Thus, 
the change of variables ( 13. 6 j) is possible for any B in evaluating (-^trS™) to give the result 

^tr5 n ^ =0 {n > 1). (3.23) 

When W'{(j)) is linear, ^tr5") has the same expression as the RHS of (I3.20p except the 
integrand multiplied by j^trB n . It leads to a finite result of the ^-integration for n > 1, 
and ( I3.23P is also obtained. 

Furthermore, it can be similarly shown that the unnormalized expectation values of 
multi-trace operators Yli=i ^tri?™ 1 (m, ■ • ■ , n-k > 1) vanish: 

k . \' 

JJ trB n. ,\ =(J _ (324) 
,i=l / 

3.1.2 Localization to W'((f) = 0, and localization versus Vandermonde 
Since (13.241) means 

°°i / i \ n i / 1 \ n\ l 



SM-^) ((> B2 )> =<1> ' = ^° (3 ' 25) 



for an arbitrary parameter u, we may compute 7Vtr ^ 2 B ) \ to evaluate the partition 

function Z a=0 . It is independent of the value of u, so u can be chosen to a convenient 
value to make the evaluation easier. 
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Taking u > and integrating B first, we obtain 



Z a =o 



■1) 



iV 2 



^V^^e-^t^'W]. (3.26) 



Then, let us consider the u — > limit. Localization to W'(0) = takes place because 

N 2 

e -mr[±W(*r) = (27r) ^ Y[ S (W'(<p) a ). (3.27) 



lim I — 



0=1 



It is important to recognize that W (0) a = for all a implies localization to a continuous 
space. Namely, if this condition is met, W'(U^<j)U) a = for v £7 6 SU(N). Thus the 
original SU(N) gauge symmetry in the matrix model makes the localization continuous 
in nature. This is characteristic of SUSY matrix models. 

The observation above suggests that in order to localize the path integral to discrete 
points, we should switch to a description in terms of gauge invariant quantities. This 
motivates us to change the expression of to its eigenvalues and SU (N) angles as 



U 



Ai 



V 



\ 



A 



U\ U e SU(N). 



(3.28) 



nJ 



This leads to an interesting situation, which is peculiar to SUSY matrix models and is 
not seen in the discretized SUSY quantum mechanics. For a polynomial W'{4>) given by 
(I2.24p . the partition function (I3.26P becomes 



J a=0 



k-l 



E^E 

.fc=l 1=0 



Lk-e-i 



(3.29) 



after the Grassmann integrals. Note that the N 2 x N 2 matrix Y^=i 9k 



the eigenvalues Y%=i 9k T,i=o A i A j 
can be expressed as 

p fc-i 



has 

[i, j = 1, ■ ■ ■ , N). Thus, the fermion determinant 



det 



E^E' 

.fc=i e=o 



ik-i-i 



fe-i 



AT 

n 

n^wl n 



,i=l 

The measure given in ( 13. 4p can be also recast to 

N 



i>j 



W'(Xi) - W'(Xj 
Aj — Xj 



. (3.30) 



d N2 (j) = C N (j[dX,^ A(X) 2 dU, 

i=l 



(3.31) 



17 



where A (A) = ni>,(Ai — A?) * s the Vandermonde determinant, and dll is the SU(N) 
Haar measure normalized by J dU = 1. Cn is a numerical factor depending only on N 
determined by 

N 

(n^)A(A) 2 e-^-K (3.32) 

i=l 

Plugging these into (13.29p . we obtain 

J i=i \i=i J U>i 

N_ 

xfi)' e -^iW_ (3 .3 3) 

In this expression, the factor in the second line forces eigenvalues to be localized at the 
critical points of the superpotential as u — > 0, while the last factor in the first line, 
which is proportional to the square of the Vandermonde determinant of W'(\i), gives 
repulsive force among eigenvalues which prevents them from collapsing to the critical 
points. The dynamics of eigenvalues is thus determined by balance of the attractive force 
to the critical points originating from the localization and the repulsive force from the 
Vandermonde determinant. This kind of dynamics has not been seen in the discretized 
SUSY quantum mechanics discussed in the previous section. 

To proceed with the analysis, let us consider the situation of each eigenvalue A, fluc- 
tuating around the critical point 4>c,f 

\ = (f>ci + Vu\i (i = 1, • • • , N), (3.34) 

where Aj is a fluctuation, and C) i, • ■ ■ , (p c ^ are allowed to coincide with each other. Then, 
the partition function (13. 33p takes the form 




2 



xe-^ESU + {y/u). (3.35) 

Although only the Gaussian factors become relevant as u — > similarly to the N = 1 case 
(I2.43p . there remain N(N — l)-point vertices originating from the Vandermonde determi- 
nant of W'(Aj) which yield a specific effect of SUSY matrix models. Before computing 
(I3.35P for a general case, let us consider the following two simple cases. 

Gaussian case In the case of W'{4>) = gi4>, where the corresponding scalar potential 
\W ((f)) 2 is Gaussian, the critical point is only the origin: C)1 = • • ■ = c Ar = 0. Then, 
( I3.35P is reduced to 

r N 9 
Z a=0 = C N (Jjdty gf J] (Ai " ^) e-^-^S (3.36) 

i=l i>j 
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where no 0(y/u) term appears since W'(<f>) is linear. By using (I3.32p we obtain the result 
nothing but eq. (B.3) in [5J: 

Z a=0 = (sga( 9l )) N2 = (sgn( 9l )) N . (3.37) 

Double-well case For W'(<f>) = g(<p 2 — fj 2 ) with /x > 0, which gives a scalar potential 
of double-well shape, each of 4>c,i is equal to fi or —fi. Let us consider the case that the 
first v + N eigenvalues are around ft and the remaining u^N around —ft: 



N — A 4 , 4>c,u+N+l — ■ ■ ■ — Cj at — —fi, (3.38) 



where the filling fractions v+,v_ satisfy v + + V- = 1. Let Zr v+iV _\ be a contribution to 
the partition function Z a=0 from small fluctuations around (13.38[) . Then, 



Za =°- zZ ( u . n\uu m z ^-y (3 - 39) 

Since 



nV:0 (, + iV)!(,_iV)!^-)- 



W"{<j> c , i ) = 2gn (z = l,... ,v+N), 

W'ted) = -2gft (i — u + N + 1, - • • ,N), (3.40) 

we have 



, N 



i=l 

. i / ■ 2 - / - - \ - 

X 



n n 

v+N>i>j>\ N>i>j>u + N+1 

x II (a« + A 3 ) 2 + 0(^). (3.41) 

N>i>v+N+1, v+N>j>l 

Note that flipping the sign A, — > — Aj for i = z/ + iV + 1, - - • , iV makes the factors in the 
second and third lines combined to the square of the single Vandermonde determinant 
A(A) 2 . Thus, Z{y + ,vJ) can be expressed by the partition function of the Gaussian SUSY 
matrix model with g± = 2gfi : 

Z {v+ , v _) = {-lY- N {2gf,) N2 C N /(n^A^e-^xi^^+OCVS) 

i=l 

= (-iy- N (sgn(2gf,)) N + 0(^). (3.42) 

Here, let Z G ,v± be the partition functions of the Gaussian SUSY matrix models with 
the matrix size v±N x v±N describing contributions from Gaussian fluctuations around 
the minima <ft = ±/i, respectively. Since 

Z G>V+ = (sgn(2^/i))^ 7V , Z Gv _ = (sgn(-2g fi)Y~ N , (3.43) 
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we can show 



J G,v + ^G,u- 



(3.44) 



in the limit u — > 0. It holds for arbitrary N, and leads to the statement (4.17) in 
the previous paper [6] in the large- N limit. Note that the integrand in the first line of 
(I3.42p cannot be factorized into the products of two functions - one is a function of \ 
(i — 1, • • • , v + N) and the other of \ {i = v + N + 1, • • ■ , N) - due to the Vandermonde 
determinant. Nevertheless, the factorization (I3.44p takes place at the level of the partition 
function. It is interesting to get more insight about the factorization, which will be useful 
to make deeper our understanding on the structure of SUSY matrix models. 
Finally, we find that the total partition function (I3.39P vanishes: 



J a=0 



N 

E 

u-N=0 



(u + N)l(y_N)\ 



sgn(2^)) w (1 + (-1))" = 



(3.45) 



which is expected from the spontaneous SUSY breaking in the case of double-well scalar 
potentials at finite N. 



General case Now, let us evaluate (I3.35P for a general superpotential. We change the 
integration variables as 

A, = l . r yu (3.46) 



then the integration of Aj becomes d\i 
(I3.35P is computed to be 



i 



A' 



Z a =o 



En 



W"( 



i=l 
N 



W"( 



c 



N 



N 



Yl dy '> 



J^oo dyi ■ ■ ■ ■ In the limit u — > 0, 



i=i 



i=l 



-i N 



^ sgn(^"(0 c )) 

<t> c :W'(<j> c )=0 



(3.47) 



Note that the last factor in the first line of (13.471) is nothing but the partition function 
of the Gaussian case with g% — 1. The last line of (I3.47P tells that the total partition 
function is given by the iV-th power of the N — 1 case f)2.40p 0. 

Furthermore, we consider a case that the superpotential W(<j)) has K nondegenerate 
critical points ai, • • • , ax- Namely, W (a/) = and W" (a/) ^ for each / = 1, • • • , K. 



9 Since X^ c . w'{<t> c )=o s S n (W"((f> c )) is equal to the mapping degree (|2.36p from (|2.41l) . Z a=0 is also 
expressed as the iV 2 -th power of the N = 1 case. It is analogous to (|A.27I) in the discretized SUSY 
quantum mechanics with the identification TV 2 = T. 
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The scalar potential \W\cj)) 2 has K minima at (f) = a\, ■ ■ ■ , ax- When N eigenvalues are 
fluctuating around the minima, we focus on the situation that 
Aj (i — 1, • • • , i>\ N) are around (p = a\ 
A Wi jv+j ( i — 1) ■ • • , V2N) are around = ai 

\ UlN+ ... +VK _ lN+i (i — 1, • • • , z^A/") are around 4> = a K , 

where ^1, •• • , Pk are filling fractions satisfying £)j=i ^ = 1- Let Z^,...^) be a contri- 
bution to the total partition function Z a= o from the above configuration. Then, 

" An 

Za =°- ^ ^Ny.-.-ivKNy. 2 ^'-^- (3 - 48) 

(The sum is taken under the constraint Y^f=i 1/ i = !•) Since is equal to the 

second line of (I3.47P with <fi C)i fixed as 

4>c,l = ■ ■ ■ = (pcviN — °1> 
<Pc,viN+\ = • • • = c ,!/iA r +i/ 2 A r = a 2; 

0c,i/iAf+---+y A '_iAf+l = • ■ ■ = Cj at = ax, (3.49) 

we obtain the generalization of the double-well case (13.441) : 

K 

Z^,-,u K ) =l[Z Gui , Z G>VI = (sgn {W" \a I ))) VlN . (3.50) 
1=1 

Zg, Vi can be interpreted as the partition function of the Gaussian SUSY matrix model with 
the matrix size viN x uiN describing contributions from Gaussian fluctuations around 
= aj. 



3.2 case 

In the presence of the external field a, let us consider ~E a {B) in (13. 9p with the action 
(I3.13P obtained after the change of variables (13.61) . Using the explicit form of the measure 
(IXTUj) and (IXTTj) . we obtain 

E a (B) = (e ia - 1) ( ~^ 2-1 Jd N2 4> (d N2 ^d N2 ^j e -Ntri±B* + iBW>w^QW>w] 

xNtv(B^) Ntr(B^) e -(^-i)^(^) ; ( 351 ) 

which is valid for B ^ 0. Although we can proceed the computation further, it is more 
convenient to invoke another method based on the Nicolai mapping we will present in the 
next section. 
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4 [e %a — l)-expansion and Nicolai mapping 

In the previous section, we tried to compute the partition function Z a in the presence of 
the external field. We have seen that the change of variables is useful to localize the path 
integral, but in the q ^ case the external field makes the localization incomplete and 
the explicit computation somewhat cumbersome. In this section, we instead compute Z a 
in an expansion with respect to (e lct — 1). For the purpose of examining the spontaneous 
SUSY breaking, we are interested in behavior of Z a in the a — > limit. Thus it is expected 
that it will be often sufficient to compute Z a in the leading order of the (e ia — l)-expansion 
for our purpose. 



4.1 Finite N 

Performing the integration over fermions and the auxiliary field B in (I3.3P with W'(<f>) in 
fl22U, we nave 

(e ia - 1)1 ® 1 + 0* Z)^ ® ^ M ) e- N *% w 'W\ (4.1) 
fe=i e=o / 

Hereafter, let us expand this with respect to {e %OL — 1) as 

N 2 

Z a = Y,(e ia -l) k Z a , k , (4.2) 



k=0 



and derive a formula in the leading order of this expansion. The change of variable <ft as 
(EEZBD recasts gj) to 



N N / p k-l 



Z a = C N [(jjd^AW J] U^-l + ^^E^-Me-^^W (4.3) 

i=l i,j=l \ k=l £=0 J 

after the SU(N) angles are integrated out. Crucial observation is that we can apply the 
Nicolai mapping for each % even in the presence of the external field 

A, = (e JQ -l)A 4 + H/ , (A 4 ), (4.4) 

in terms of which the partition function is basically expressed as an unnormalized expec- 
tation value of the Gaussian matrix model 

r N 

Z a = C N (n^) II( A * - A,)V^ A ?e-^(-^+^ 2 ), (4.5) 
J i=i i>j 

where A = e ia — 1. However, there is an important difference from the Gaussian matrix 
model, which originates from the fact that the Nicolai mapping f)4.4p is not one to one. As 
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a consequence, Aj has several branches as a function of Aj and it has a different expression 
according to each of the branches. Therefore, since the last factor of ( 14. 5 p contains Aj(Aj), 
we have to take account of the branches and divide the integration region of Aj accordingly. 
Nevertheless, we can derive a rather simple formula at least in the leading order of the 
expansion in terms of A owing to the Nicolai mapping (I4.4p . In the following, let us 
concentrate on the cases where 

Aj — > oo as Aj — >■ ±00, or Aj — > —00 as Aj —> ±00, (4.6) 

i.e. the leading order of W'(4>) is even. In such cases, we can expect spontaneous SUSY 
breaking, in which the leading nontrivial expansion coefficient is relevant since the zeroth 
order partition function vanishes: Z a=Q = Z a $ = 0. Namely, in the expansion of the last 
factor in (14. 5 p 

e -;vEf =1 (-AA,A i+ i^) =1-N^2 + \A 2 \ 2 ^j + • • • , (4.7) 

the first term "1" does not contribute to Z a . It can be understood from the fact that 
it does not depend on the branches and thus the Nicolai mapping becomes trivial, i.e. 
The mapping degree is zero. Notice that the second term also gives a vanishing effect. 
For each i, we have the unnormalized expectation value of N (AAjAj — ^ A 2 A?) , where 
the Aj-integrals (j 7^ i) are independent of the branches leading to the trivial Nicolai 
mapping. Thus, in order to get a nonvanishing result, we need a branch- dependent piece 
in the integrand for any Aj. This immediately shows that in the expansion (14. 2p . Z a ^ = 
for k = 0, . ■ ■ ■ , N — 1 and that the first possibly nonvanishing contribution starts from 
0(A N ) as 



N N 



Z a , N = C N N N J (jjdK) n(Ai- Aj)>e-»^W I1(A 4 A 

i=l i>j i=l 



A=0 



Note that the A = (e ia — Independence of the integrand comes also from Aj as a function 
of Aj through ( 14. 4p . Although the integration over Aj above should be divided into the 
branches, if we change the integration variables so that we will recover the original Aj 
with A = (which we call Xj) by 

Aj = Wfa), (4.9) 

then by construction the integration of Xj is standard and runs from —00 to 00. Therefore, 
we arrive at 

Z a , N = C N N N (T[dxA]\{W''{x i )W'{x i )x t )\[{W ! {x t )-W ! {x j )) 2 

J -co \_1 ' „-_1 „-^„- 



i=l i=l i>j 

xe-tfESa^W (4.10) 
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which does not vanish in general. For example, taking W(<f>) = g(<p 2 — p 2 ) we have for 
N = 2 

Io 5 V/n 



Z a , 2 = Wg 2 C 2 ll 



(4.11) 



where 



dXX n e- 9 (x -^ (n= 0,2,4, ••-)■ 



(4.12) 



In fact, when g = 1, p? = 1 (double- well scalar potential case) we find 

2 



Io = 1.97373, 

5 Wo 



-0.165492 ^ 0, 



(4.13) 



hence Z ai 2 actually does not vanish. In the case of the discretized SUSY quantum mechan- 
ics, we have seen in flA.36j) that the expansion of Z a with respect to (e ta — 1) terminates 
at the linear order for any T. Thus, the nontrivial 0(A N ) contribution of higher order 
can be regarded as a specific feature of SUSY matrix models. 

We stress here that, although we have expanded the partition function in terms of 
(e ia — 1) and (14.101) is the leading order one, it is an exact result of the partition function 
for any finite N and any polynomial W'(<p) of even degree in the presence of the external 
field. Thus, it provides a firm ground for discussion of spontaneous SUSY breaking in 
various settings. 



4.2 Large-iV 



As an application of (14.101) . let us discuss SUSY breaking/restoration in the large- limit 
of our SUSY matrix models. From (I4.10p . introducing the eigenvalue density 



1 N 

i=i 



X ; 



the leading 0(A ) part of Z a is rewritten as 



r N 

Z a ,N = C N N N (n^) exp(-iV 2 F) 
J i=i 



with 



N 



dxdyp(x)p(y) log \W'(x) - W\y)\ + / dxp(x)-W (x) 2 



dxp(x)\og(W"(x)W'(x)x). 



(4.14) 



(4.15) 



(4.16) 
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In the large-iV limit, p(x) is given as a solution to the saddle point equation obtained 
from 0(1) part of F as 



provided that there exists an 0(1) solution of this equation. Plugging a solution po(x) 
into F in (14.161) . we get Z a in the large- N limit in the leading order of (e ia — l)-expansion 

as 

Z a , N ^N N exp(-N 2 F ), 

Fo = -J dxdyp (x)p (y) log \W'(x) - W'(y)\ + J dxp (x)^W'(x) 2 - ^\ogC N , 

(4.18) 

where Cn is a factor dependent only on N which arises in replacing the integration over 
(f) by the saddle point of its eigenvalue density, thus including Cn- From consideration of 
the Gaussian matrix model, Cat is calculated in [6] as 



C N = exp 



-N 2 + 0(N°) 
4 



(4.19) 



and is expected to be independent of the form of superpotential. In (I4.18P we notice 
that, if we include 0(1/N) part of F (the last term in (I4.16P ) in deriving the saddle 
point equation, the solution will receive an 0(1/ N) correction as p(x) = po(x) + -^pi(x). 
However, when we substitute this into ( I4.16|) . pi(x) will contribute to F only by the order 
0(1/ N 2 ), because 0(1/ N) corrections to F under po(x) —> po(x) + -hpi(x) vanish as a 
result of the saddle point equation at the leading order (I4.17P satisfied by po(x). 
On the other hand, if we set a = at the level of f)4.3p . we have 



r N N 

• , = / (11^) \[W"(\) WiWW) - W'(\ 3 )) 2 e- N ^^ w '^\ (4.20) 

8=1 i=l i>j 



from which we obtain exactly the same saddle point equation as (I4.17P at iV = oo. Namely, 
making the expansion with respect to (e ia — 1) affects only the subleading part of F in the 
1/iV-expansion as one can see by comparing (14. 10p and (I4.20p . It is also the same as the 
saddle point equation (3.15) in the previous paper [6]. Thus, various large-iV solutions 
derived in section 4.1 in |6j, which restore SUSY, can be reproduced from Z a ^ or Z a in 
the large- N limit followed by the a — > limit, in spite of the SUSY breaking at any finite 
N (Z a=0 = 0). Let us see it explicitly for the free energy. When N is large but finite, the 
twisted partition function will take the form 

Za ~ (e ia - 1) N N N c e - N2F °- NF ' , (4.21) 
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where Ft is 0(1/ N) contribution of F given as 

F x = - J dxpo(x) log(W"(x)W'(x)x), (4.22) 

and the coefficient c comes from 0(1/ N 2 ) contribution of F. The free energy, which 
corresponds to the quantity —-^2 log|Z a | leads to F Q in the large- N limit followed by 
a — > 0. Notice that, although the effect of the twist (e ia — 1) N is of the subleading order 
at large N, it plays a crucial role to obtain the large- N free energy F . (If a was sent to 
zero before the large- N limit, we would have the vanishing partition function and could 
not find the large- N free energy F .) 



4.3 Example: SUSY matrix model with double-well potential 

For illustration of results in the previous subsection, let us consider the SUSY matrix 
model with 

W'((j>) = <p 2 -fi 2 (fi 2 G R). (4.23) 
In this case the saddle point equation ( 14.1 7p becomes 

p(y) , f p(y) _ 3 .2 



x — y J x + y 



x 6 - [fx. (4.24) 



In section 4.1 of [6], we have obtained an asymmetric one-cut solution where the eigenvalue 
density has a single support [a, b] with b > a > and also a two-cut solution with a 
symmetric support [—6, —a] U [a, b] 0. Here, a 2 = — 2 + /j 2 , b 2 = 2 + /i 2 , thus they 
are valid for fi 2 > 2. In N — 1 case, it is well known that the SUSY is spontaneously 
broken for (14.231) . but in the matrix model case we have shown in [6] that the SUSY 
is restored in the large-iV limit for both solutions. In particular, the free energies for 
both solutions are shown to vanish and therefore they coincide with the value of the 
free energy of the Gaussian matrix model. It is also proven that the expectation values 
(-^tri? n ) (n = 1, 2, • • • ) are all nil. Here it is worth pointing out that the principal value 
in the saddle point equation, in particular in the second term in ( I4.24p plays a crucial 
role in the existence of the two-cut solution. In this subsection we investigate a one-cut 
solution with a symmetric support [— c, c] which has not been discussed in [B]. 



10 Note that /i 2 here corresponds to — /j 2 in section 4.1 of [6 r Interestingly, the eigenvalue distribution 
of the two-cut solution is not Z2 symmetric in general. In fact, 



= ftV&EE W^ ^< X <^ (4.25) 
FK ' \^MV(z 2 -a 2 )(6 2 -x 2 ) (-b<x<-a) V ' 

is the explicit form of the solution with the filling fraction {v +1 ^_), which includes the asymmetric one-cut 
solution as a special case {v + , v ) = (1,0). 
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At first sight, it seems strange that there exists such a solution because the fermion 
determinant in the partition function (14. 3 p looks 

N N 

Y[ ( e ia - 1 + Ai + Xj) = J] {e ia - 1 + 2A,) J] (e ia - 1 + A, + A,) 2 , (4.26) 

i,j=l i=l i>j 

the first factor of which makes Aj apart from the origin in the a — > limit. However, as 
we will see below, the symmetric one-cut solution exists owing to the large- N limit and 
we will confirm its validity by checking finiteness of the free energy for our solution. 

In order to solve (14.24p for p(x) with a symmetric support [— c, c], let us consider a 
complex function 

G(z) ee f dy^, (4.27) 



and further define as in H3 



z-y 



G_{z)= l -{G{z)-G{-z)), (4.28) 

then G_(z) has following properties: 

1. G-(z) is odd, analytic in z G C except the cut [— c, c]. 

2. G-(x) G R for x G R and x ^ [— c, c]. 

3. GL(z) ->■ i + 0(4f) as z ->■ c». 

4. G_(x ± iO) = |(x 2 — yU 2 )x =F iirp(x) for x G [— c, c]. 
They lead us to deduce 

G-{z) = \{z 2 - p 2 )z -\U~» 2 + y) (4.29) 



with 

from which we find that 



2tt 



c 2 (3c 2 -4/z 2 ) = 16, (4.30) 



fx 2 - p 2 + ^0 Vc 2 - x 2 , x G [-c, c\. (4.31) 



Po^) > tells us that this solution is valid for p 2 < 2, which is indeed the complement 
of the region of p 2 where both the two-cut solution and the asymmetric one-cut solution 
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obtained in [6] exist. Given po(x), it is straightforward to calculate the free energy ( I4.18P 
as 

F (x) = \x 2 - -^-x 4 - — (x 3 + 30x) Vx 2 + 12 - log(x + Vx 2 + 12) + log 6, (4.32) 
3 216 216 

where x = p 2 . In contrast to this, it is observed in [6] that for p 2 > 2 the free energy 
calculated from the asymmetric one-cut solution or the two-cut solution is independent 
of p 2 and vanishes, reflecting the restoration of the SUSY. It is easy to see that F (x) > 
for x < 2, and the expectation value of j^trB is computed to be 



c 



i 



,,2 



— (c — x) — X 
16 V ; 



(4.33) 



which is nonzero for x < 2. These are strong evidence suggesting the spontaneous SUSY 
breaking. Also, the x- derivatives of the free energy, 

lim F (x) = lim Hm = , Hm = --, (4.34) 

x~>2-0 dx x^2-0 dx 2 x-^2-0 dx 3 2 

show that the transition between the SUSY phase (x > 2) and the SUSY broken phase 
(x < 2) is of the third order. 

As commented in ( I4.26p . if we take a look at the 0(1/N) contribution of F given in 
(I4.22p . naively it seems strange that we have a nonzero saddle point eigenvalue density 
around the origin, since xW"(x) = 2x 2 ~ there and the integrand of f)4.22p diverges. 
Furthermore, it is also curious that we have eigenvalues in general distributed around 
zeros of W'(x), because the integrand again diverges there which would mean that the 
partition function vanishes. However notice that, even if the integrand looks divergent, it 
is just logarithmic and its integral itself is finite due to a contribution from the measure. 
Because Pp(x) is finite, the relevant integral over the vicinity of the singularities for the 
real part [Jj of F\ is at most J Q e dx log x which clearly converges at the origin. (The 
logarithmic singularity is integrable.) Therefore, 0(1/N) part of F does not diverge 
owing to the large- N limit. 

As an example, in the double-well case W'(<p) = 4> 2 — fi 2 with fi 2 > 2, let us consider 
the two solutions obtained in [B]. The asymmetric one-cut solution is given by 



Po(x) 



-y/(x 2 -a 2 ){b 2 -x 2 ), a 2 = -2 + p 2 , b 2 = 2 + fi 2 , (4.35) 



7T 



and the two-cut solution with the filling fraction (z/ + , z/_) is (14.251) . Evaluating (I4.22p at 
each of the solutions, we find the same result of the real part of F\ for both solutions : 

Re F 1 = - log (V + vV - 4) + 1 - i/i 2 (V - vV " 4) . (4.36) 



11 The imaginary part is irrelevant in the analysis, because it just contributes to the overall sign of 
the partition function as sgn ^FJi^i {W'(xi)W'(xi)xi) \ in evaluating the partition function at a single 
large- N solution. 
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It is finite as understood from the above reasoning, and interestingly it is not dependent 
on u±. The finiteness supports the validity of our large- N solutions. 

Finally, we make a comment on one of interesting aspects of our result that the exis- 
tence of a leading nontrivial contribution at O ((e* Q — 1)^) given in (14.16P suggests SUSY 
breaking for finite N, but that in the double-well case (I4.23P with fj, 2 > 2 it leads to the 
supersymmetric solutions to (I4.17P in the large-iV limit with vanishing free energy, as 
discussed below (I4.20p . 

5 Summary and discussion 

In this paper, firstly we discussed localization in discretized SUSY quantum mechanics 
without the external field a by changing integration variables. It makes it clear that the 
path integral is localized at the auxiliary field B = 0, which in turn implies the standard 
localization at the critical points of superpotential. Furthermore, it was investigated in 
detail how Eq(B), the integrand of the partition function with respect to B, behaves as 
B ~ 0, and clarified whether the change of variables is applicable or not. Similar argu- 
ments were presented also for a ^ case. This gives a different approach to localization 
from a deformation by Q-exact terms. It is worth pointing out that the change of variables 
can be applied even to systems where SUSY is (explicitly) broken, while a deformation 
by Q-exact terms cannot be straightforwardly. Thus, the former is useful to investigate 
localization in systems where the external field explicitly breaking SUSY is turned on. 
We also stress that we provided a firm formulation of change of variables for localization 
(without issues mentioned in footnote [3]) in the path integral which is useful in discussion 
of spontaneous SUSY breaking. As emphasized in the introduction, such a formulation 
is indispensable because nonperturbative formulations of string theory in terms of matrix 
models are defined by the path integral. 

Secondly, we explained localization in SUSY matrix models without the external field. 
The formula of the partition function was obtained, which is given by the iV-th power 
of the localization formula in the N — 1 case (N is the rank of matrix variables). It 
can be regarded as a matrix-model generalization of the ordinary localization formula. In 
terms of eigenvalues, localization attracts them to the critical points of superpotential, 
while the square of the Vandermonde determinant originating from the measure factor 
gives repulsive force among them. Thus, the dynamics of the eigenvalues is governed 
by balance of the attractive force from the localization and the repulsive force from the 
Vandermonde determinant. It is a new feature specific to SUSY matrix models, not seen 
in the discretized SUSY quantum mechanics. For a general superpotential which has K 
critical points, contribution to the partition function from ujN eigenvalues fluctuating 
around the J-th critical point (J = 1, • • • ,K), denoted by Z( Vlt ... tl/K ), was shown to be 
equal to the products of the partition functions of the Gaussian SUSY matrix models 
Zg, Ui ■ ■ ■ Zg,» k - Here, Zg, Ui is the partition function of the Gaussian SUSY matrix model 
with ujN the rank of matrix variables, which describes Gaussian fluctuations around 
the 7-th critical point. In the double-well case, it leads to the claim of the previous 
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paper [6J. It is interesting to investigate whether such a factorization occurs also for 
various expectation values. 

Thirdly, as mentioned in the above, the argument of the change of variables leading 
to localization can be applied to a ^ case. Then, we found that a-dependent terms 
in the action explicitly break SUSY and makes localization incomplete. Instead of it, 
the Nicolai mapping, which is also applicable to the a ^ case, is more convenient for 
actual calculation in SUSY matrix models. In the case that the supersymmetric partition 
function (the partition function with a = 0) vanishes, we obtained an exact result of a 
leading nontrivial contribution to the partition function with a ^ in the expansion of 
(e ia — 1) for finite N. It will play a crucial role to compute various correlators when SUSY 
is spontaneously broken. Large- N solutions for the double- well case W'(4>) = (ft 2 — fi 2 
were derived, and it was found that there is a phase transition between the SUSY phase 
corresponding to /i 2 > 2 and the SUSY broken phase to fi 2 < 2. It was shown to be of 
the third order. 

For future directions, this kind of argument can be expected to be useful to investigate 
localization in various lattice models for supersymmetric field theories which realize some 
SUSYs on the lattice B 

Also, it will be interesting to investigate localization in models constructed in ref . , 
which couple a supersymmetric quantum field theory to a certain large- N matrix model 
and cause spontaneous SUSY breaking at large N. 

Finally, we hope that similar analysis for super Yang-Mills matrix models [21 [31 H], 
which have been proposed as nonperturbative definitions of superstring/M theories, will 
shed light on new aspects of spontaneous SUSY breaking in superstring/M theories. To 
carry out the analysis, the method of the Gaussian expansion or improved perturbation 
theory would be useful [261 El EH EH [29]. 
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A Localization in discretized SUSY quantum mechan- 
ics with general T 



In this appendix, we generalize the argument for T = 1 presented in sections 12.21 and 12.31 
to T > 2 case. 

A.l Localization for general T 

The action (12. 9p with the periodic boundary condition (a = 0), denoted by So, is invariant 
under the Q-SUSY: 

Q0(t)=^(t), Q^{t) = Q, Q$(t) = -iB{t), QB(t) = 0. (A.l) 

We consider the field redefinition 

0(t) = 4>(t) + e^(t), m = kt) - ieB(t), (A.2) 

where ip(t) is chosen to be orthogonal to —ieB(t) as 

T 

Y^B(t)$(t) = 0. (A.3) 
t=i 

Namely, the degrees of freedom of ip{t) are carried by both e and ip(t). Differently from 
the T = 1 case, e cannot parametrize the whole functional space of ip(t), and represents 
merely a single degree of freedom of ^(t) parallel to B(t). 

Similarly to the T = 1 case, since flA.2j) is recast as a shift by the Q transformation: 

cj>{t) = fa) + eQ4>(t), m = $(t) + eQ$(t), (A.4) 
the action So is shown to be independent of e: 

So{B,(j)^^) = S {B + eQB, 4> + eQ4>, + eQ^fj, $ + eQi/j) 
= S (B,<p,ij>,ij>) 
T r l 



^B(t) 2 + iB(t) \ <j>{t + 1) - 0(t) + W'm)) 



+ kt) fa + 1) - m + w"m))m}} , (a.5) 

from the Q-SUSY invariance of So- 



31 



When we write the partition function as 

T 

Z = I f[dB(t)Z (B), (A.6) 



T P T 

So 



E (B) = j n^(t) dm d^t)) e 

= (-l) T(T - 1)/2 /n#w(^J [lKdcf>(t)dm)e- S °, (A.7) 

we can regard B(t) and ip(t) as external fields in the integral J ^n^=i d<p(t) d^)(€) \ e~ s ° . 
This point of view makes easier to derive the Jacobian of the path-integral measure 
associated with (1A.2j) . For —ieB(t) in (IA.2j) . we decompose it as 

- ieB(t) = -iAf B e x -r^-B(t), (A.8) 



with Mb = = Jj2t=i B(t) 2 - Since the second factor j^B(t) can be regarded as "a 
normalized wave function" , i.e. a unit vector in the functional space of B(t), the remaining 
—iNb£ is identified with an integration variable. Thus, the measure of e associated with 
flA.2j) is given by <i(— i/Vse) = jj^de. Also for ip(t), expressing the constraint flA.3j) as a 
delta-function, the measure is explicitly defined by0 

T / T N 

T „ T 



= (-if- 1 (n £ (a.9) 

Hence, we obtain the measure for xjj as0 

T 

J[dm = j^de (d^Y (A.10) 
After the change of variables flA.21) . Ho(-B) becomes 



T 

x / / * e-Sr=^ = (*){^(*+ 1 )-^(*)+^"(^(*))V'(*)}. (All) 

t=i ^ 



13 The sign factor (— 1) T 1 can be determined so that the RHS becomes Yit-2 when B(t) vanishes 
except B(l). 

14 For T = 1, fd?^ is reduced to -^p, so we have = ^3- S(l)c2e = -g^y de which reproduces the 
Jacobian in (I2.20j) . 
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It is clear that S (-B) vanishes due to the trivial e- integral as long as Mb 7^ 0. Thus, the 
path integral of the partition function is localized to Mb = 0, i.e. B(l) = ■ ■ ■ = B(T) = 0. 

On the other hand, as in (12.281) of the T = 1 case, the change of variables flA.2j) does 
not work when Mb ~ 0. In order to see how E (B) becomes singular as Mb — > 0, let 
us consider flA.llj) in the region Mb < £• We express (-B(l), ■ ■ • ,B(T)) e R T by polar 
coordinates with the radial direction Mb an d the angular directions specified by the unit 
vector 

Q B (t) = ^-B(t). (A.12) 
For W'(4>) given by (I2.24p with p > 2, we rescale 4>(t) and ip(t) as 

^ = (aTb) P ^ = U ^ ^ (t) ' (A ' 13) 

which correspondingly changes the measure as 

= f^-j p n #'(*). {<$) = [-j^) P W • ( A - 14 ) 

The rescaling is convenient to see the A/s-dependence of S (-B). Then, the integrands of 
the ^-integral and the Grassmann integral become 



-iT:J=iB(t){kt+^-4>(t)+W(4>(,t))} - e -iELi^(*)9 P ^'(*) p [1 + 0(e 1/p )] , 

-Yj^kt){^{t+i)-m+w"m))m} _ e -Ef=i^' (^vw [i + o^ 1 /?)] 

(A.15) 



for Mb < e, respectively. Plugging the above results, S (-B) can be expressed as 



7M 

*ELi^bW9p'/ > '(*) p 



T 

X 



n#(*) / e-^^'W^ [l + 0(e x /P)] . (A.16) 

We thus find that S (-B) becomes singular as M B p for Mb ~ 0. 

More precisely, if we define = J Afg<£ Y\J =1 dB(t)E (B) as in (12.221) . the factor 

M^ 1 coming from the measure Yl[ = i dB(t) = M B ^ l dMsdVtB expressed in the polar 
coordinates makes somewhat milder the singularity at Mb — in H (-B). However, it is 
not sufficient to achieve convergence because 

f « T_i (£f - [ m ° (£f* x i 1 + o<£2) ] = ~- < a - i7) 
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Hence, in Z^°\ the trivial e-integral vanishes while the TV^-integral diverges. We see that 
the change of variables (IA.2I) is not appropriate for , because it leads to an indefinite 
expression: = oo x 0. (In p = 1 case, we can consider W'((j)) = because the 
constant term go can be absorbed by a shift of <p. The ^-integrals in flA.111) yield 



T 

, \f / fJd0(t) e - i E"=i fl W{^+ 1 )-^*)+ w 'W*))} 
(2tt) J t=1 

f[S ((1 - - B(« - 1)) = ,, , T| t[sm)) (A.18) 

*=i I ^ 9i) \ t=1 

with 5(0) = B{T). Due to the delta-functions Y^ =1 5(B{t)) and in 1KH1) . the 

5-integrals in Zq '' become singular leading to an indefinite form of Zq ^: oo x 0. ) As in 
the T = 1 case, the indefinite form of Z^ obtained after the change of variables (1A.2I) 
implies that Zq ^ possibly takes a nontrivial value. 

Unnormalized expectation values For the unnormalized expectation values of B(t) n 
(n > 1): 

(B (t) n y = I (II rf5 (*)) B (t) n So(B), (A.19) 

by the same change of variables, the contribution from the region Mb — \ \B\\ > s clearly 
vanish from J de = 0. Hence, 



e < ELi{( 1 -9i) B W-- B ( t - 1 )}^( i ) 



/ i\T(T+l)/2 



(27T) 



B||<e 



T T 

t=i b j t=i 



x j de (#) e- 30 ^'^. (A.20) 

If is given by f !2.24p with p > 2, using (jA.160 and the polar coordinates for 

(5(1),- •■ ,5(T)), we have 

= i dMsM^e-^ Yn (t) [I + O(eVP)] , (A.21) 
rffi B fi B (t) n / JJrf0'(t)e- iE ^ ! 



(_i\T(T+l)/2 



(2*0 T 

T 

X 



t J =1 f2 s (t) 9p 0'(t)P 



/ n#(^) / de(d$) e -S^^(*)MP^(*)*-V(t). (A.22) 
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Note that, since the A/s-integral is not singular at the origin for n > 1, the field re- 
definition ( ]A.2j) is always possible, differently from the case of the partition function. 
(For completeness, we show that Y n (t) definitely vanishes, i.e. factors in front of trivial 
Grassmann integrals are finite (in fact, they vanish) in appendix IA.3I ) Thus, we can show 

(B{t) n )' = (n > 1), (A.23) 

due to the trivial e-integral. In p = 1 case, the delta-functions YlJ = i 5 (B(tY) arise after 



the <fi integration in ( lA.20p . which makes the 5-integrals finite for n > 1 as 
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/ (fl^W) ^e"^ f[S(B(t)) = 5 nA . (A.24) 



It leads to flA.23j) from the trivial e-integral. 



In general, we find that the unnormalized expectation values of I?(l) ni ■ ■ ■ B(T) nT with 
n\, ■ ■ ■ , rir — 0, 1, 2, • • • and Ylt=i n t — 1 vanish: 

(B(l) ni ■■■B(T) nT )' = 0. (A.25) 
Localization to 0(t + 1) - 0(f) + W'((j)(t)) = Because 



e 2 



holds for an arbitrary parameter u from (IA.25D . the partition function can be computed 
similarly to the T = 1 case. Taking u > and integrating with respect to B, we have 

^0 = |n#(t)^^ e -^-^^ 

T 

x(-lf / (jj e-^^'W^( t+1 )^W +H/ "^)^*». (A.28) 

t=i 

In the limit u — > 0, the integration with respect to 4>(t) is manifestly localized to config- 
urations satisfying cj)(t + 1) - 0(f) + W'((j)(t)) = 0. 



15 The result (|A.24I) is obtained, if we integrate B(t') (V ^ t) before the S^-integral. Otherwise, 
we would have the vanishing result even for n = 1. We choose the order so that the result is reduced to 
([2321 when T = 1. 

16 An explicit computation of Zq is given in appendix A in |30j . where deformation invariance by 
Q-exact terms is used to obtain 

K <j>:W{<t>)=0 1 / 

It is the T-th power of the result of the T = 1 case. 
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A. 2 Localization in the presence of external field 

The action (12. 9 j) with the twisted boundary condition (12.101) . denoted by S a , can be 
written as 

S a = S + (e ia -l)^(T)^(l), (A.29) 

where So is the action with the periodic boundary condition. As mentioned below (12. 10 j) . 
<fi(T + 1) and ip{T + 1) were replaced with (f)(1) and e za i/;(l), respectively. Under the 
change of variables (I A . 2 [) which is defined for variables at t = 1, • • • , T, this becomes 

S a (B, 0, i(>, $ = 5 (S, 0, ^ $ + (e fa - 1)^(T)^(1) - *(e 4a - l)eS(T)^(l). (A.30) 

The first term is given by (IA.5j) . independent of e. Due to the last term —i[e ta — 
l)eB(T)ip(l), the e-integral in the partition function does not vanish. 
Similarly to the a = case, let us write the partition function as 



Z a = / ]JdB(t)E a (B), (A.31) 

rp rp 

E a (B) = (^j J f[ (dcf>(t) #(t) #(*)) e- s « 

= (-l) T{T - 1)/2 [ U d ^ (y) /nW*)^))^, (A.32) 

i=l \ n / J t=1 



also 



33) 



^ 0) =/ fTdS(t)S a (S), Z a = / fTtZB(t)S a (S). (A. 

■/||B||<e£=i ■'l|B||>e t=1 

The field redefinition (IA.2I) recasts H a (.B) to 

T 

x I (J[dm) (4) e - 5o(B '^^a)=o, (A.34) 
after integrating over e and -0(1) . In the process, we used 

#(1) dee i{eia - 1] ~ eB{TWl) = J #(1) de [l + t(e ia - l)eS(T)V>(l)] 

= i(e iQ - l)B(T). (A.35) 
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Note that it is valid for Mb — \\B\\ 7^ 0. In the case Mb ~ 0, we should keep the first 
term "1" in the expansion of e l ^ c " _1 ) £ - B ( T )^( 1 ) in the RHS of the first equality in flA.35[) . 
although it gives the vanishing ip(l)- and e-integrals. Integrating it over \ \B\ \ < e yields a 
singularity at the origin, so we have an indefinite form (00 x 0) which cannot be discarded 
safely It is parallel to the situation of the T = 1 case discussed in footnote |6j 

The contribution to the partition function from the integration region \ \B\ \ > e: Z a is 
in general nonvanishing. However, when W'(<j)) is linear, the ^-integrals yield YlJ=i $(B(t)) 
leading to Z a = 0. 

On the other hand, the partition function can be computed directly from ( 12. lip without 
using (lA.2p . Since the fermion determinant can be written as the sum of the determinant 
under the periodic boundary condition (n^li( — 1 + W"{4>{t))) — (— 1) T ) and the effect of 
the twist ( — (— l) T (e ta — 1)), we have 

Z a = Z - (e» - 1) (^=j J e -^ =1 {mi)-<Kt)+wmm\ (A.36) 

Since the second term is the net effect of the twist, it should be equal to the sum of 
Za^ — Zq QS) and Z a . Note again that although it vanishes in the a — > limit, it becomes 
important when the SUSY is spontaneously broken, i.e. Zq = 0. Let us elaborate on the 
former as in the T = 1 case. When W'{4 > ) ls linear (W'(<f>) = g\<f> + go), we explicitly 
obtain 

z (o) _ z (o) = _ {eia _ 1} _ ( ^ r , Z a = 0, (A.37) 

which again means that the localization takes place even in the presence of the external 
field, and that the effect of the twist on Z^ remains even in the e — > limit. It can 
be understood from S a (B) being proportional to Hq(S) similarly to the T = 1 case. In 
contrast, we show that 

Zf) - 2®> = -(e- - 1) (^-Y I (f[dB(tj) e-^»W 

'f[d<Ktj) e -*E^iB(t){«*+D-^)+w'Wt))} (A . 38) 



t=i 

vanishes as e — > when W'(4>) is a polynomial ( I2.24p with p > 2. After the rescaling 

1 

(pit) = j P <j>'(t), the ^-integrals become 



e -*ELi B (*){*(*+i)-0(t)+W(0(t))} 



" I (n<W)) e-^^t^m* [1 + ( e i/p)] . (A.39) 
J t=i 
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Then, 

zv>-z® = -(^- 1 )fe) r (f^ iHe ' KB ) y 

x [1 + 0{e 1/p )] , (A.40) 

T 

r = [<m B /"(n^))^^^ 3 ^'^, (A.4i) 

where the A/s-integral is 0(£ T ^ 1_ p 3 ), and K is also shown to be finite in appendix IA.4I 
Thus, we see that Za^ — zjf^ = 0(e T ^~p^) vanishes for p > 2 as e approaches to zero. 



Unnormalized expectation values After the change of variables (1A.2j) . the unnor- 
malized expectation values of B(t) n (n > 1): 

r T 

(BW n )l = / (II dB ^) E -(B) (A.42) 

J t=i 

are expressed as 

mrt = ( eto - 1 ) ( ~ 1 ( C)T 1)/2 jf[(dB(t)dm) ^B(tr 

T 

x f (j[ty(tj) {4) e- flo(fl '^U=°. (A.43) 



Here, since the S-integrals are not singular for n > 1, we can safely drop trivial Grassmann 
integrals, differently from the case of the partition function Z a . This is also the case for 
more general expectation values (£>(l) ni • • • B(T) nT )' a with n 1; • • • , rix > and Y^=i n t > 
1. 

A. 3 Computation of Y n (t) 
Y n (t) in f[02|) given as 



Y n (t) = y , ^ fl n B (i) w y"n#'(t)e- i ^= 1 ^( t )^'( t r x(0 / )) (A44) 



T 

X(0') ee (-l) r(T - 1)/a y"(n#(*)) (#0 e-^^'^W^CA^) 

has trivial Grassmann integrals with respect to e and one of ^' which give zero. Here, we 
show that factors in front of these integrals are finite (in fact, they vanish), which means 
that Y n (t) definitely vanishes. 
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First, after ^-integrals, we have 

T 



Xtf) = (-l)^-iV^ fapgrfitr 1 ) J dej (#') f|?W- (A.46) 
Using (#') = (-l) 7 ^ 1 (nf=i#'W) (EL^bW^'W) from the definition (TAT9]) . we 



obtain 

T 



X(0') = (-If 2 (n^W-'j / *(/ , (A.47) 

where the Grassmann integrals J de and J dip' (t)^' (t)ip' (t) trivially vanish. Thus, f lA.44j) 
has the form 



Y n (t) = j ^dO B B (t)-n(^/° C ' dtXQe^V^Wpg^t)*- 1 

x n B(t) j de #'(t)^(t$'(t)) . (A.48) 



For the (//(t)-integral, using the Nicolai mapping, we obtain 

1 /»oo 1 /»oo 

— / d<f)'{t) e -M B (t) 9p <p'{ty pgtftyr- 1 = jt — / rfx(t) e-^ww) = (j 5(n B (t)). 

(A.49) 

(j is the mapping degree of X(<f>) = g p (p p - 

sgn(^p) for p: odd 
for p: even, 

which coincides with the mapping degree of = g p p + g p -i(p p ~ 1 + ■■■ + go- 

As a result, F n (t) has a form 



Y n (t) = (J 



r 



(TXT 1 
fj5(0 B (t))J ^fi B (t) 1 #'(t)^(t)$'(t)J , (A.51) 



here we note that, since (fi B (l), • • • , f2 B (T)) is a unit vector in R T , the integration region 
for Qb is the unit (T — l)-sphere S T_1 and does not contain the origin, i.e. the support of 
Ylt=i <K^b(£))- Hence, the f^-integrals flA.51[) vanish, which shows that Y n (t) is definitely 
zero. 
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A. 4 Finiteness of Y 

In this appendix, we show that 



T 

y= fdn B f (n^'W)e-* ELi 

J J t=l 



n B (t)g p ct>'(t)f 



in flA.411) is finite for p > 2. 
From 



/g-fa^f = ^ | a ]l 

■ A cos I ( I i) ip: odd) 



for a G R, we have the bound for \Y\: 

m<(-V(i + i 

In the polar coordinates 

fis(l) = COS#i, 
Q B (2) = sin^i cos 9 2 , 



dn B f[\n B (t)\-p. 



t=i 



(A.52) 



(A.53) 



(A.54) 



n B (r-i) 
n B (T) 



sin 6*i sin 9 2 ■ ■ ■ sin Ot-2 cos 0t-i, 
sin 6*i sin 9 2 - ■ ■ sin T _ 2 sin 

with < #i, • • ■ , #r-2 < and < 9t-i < 27r, the measure is given by 

d£l B = sin T " 2 #i sin T " 3 # 2 ■ ■ • sin T _ 2 d9 1 d9 2 - ■ ■ d9 T _ 2 d6 T -i- 

Then, the f^-integrals in (lA.54j) can be expressed as 

T T-1 

dn B ]l\n B (t)\-p = 2 T J] 



*=i 



d0, fcos^H (sinfl/)^"^ 1- ^" 1 



(A.55) 
(A.56) 

. (A.57) 



For each t, since both of the powers of cos 9t and sin 9 t are greater than —1, the ^-integral 
is finite. Thus, ( 1A.57j) is finite, meaning that \Y\ is so. 



B Computation of Yj 



N 



Y N given in (15321) . 



-1 

v/2tt 



A f2 



^tr(^) y>Ve- imr(tW ' P) 



N 2 -IJ' -Ntr 



[i>'9 P ELo 1 1 



(B.l) 
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contains vanishing Grassmann integrals. In this appendix, we explicitly compute Y N to 
show that prefactors of the vanishing Grassmann integrals are finite. 

In terms of coefficients in the expansion by the basis {t a }, the measures are expressed 

as 



d 



N 2 il 



N 2 

n 

a=l 



2tt 



N 2 

d N2 iJ = i[diJ a , 



N 2 



a=l 

N 2 



d wa -y = (-i) jva - i (n^ te ) E fi ^ 

After the Grassmann integrals, we obtain 



(B.2) 



vO=l 



a=l 



Y: 



N 



{N 2 



N 2 



a=l 



xJ^^B dl \ > 



(B.3) 



where 



W) a = iVtr(t%0' p ), 

P-i 



/\a& 



5 



mr ( r^^^v^- 1 ) = ^wr- 



^=0 



(B.4) 



V{((f)') gives the Nicolai mapping to recast the (//-integrals to J (jj^=i ~7|=) e l ^"= in s y i a . 

The mapping degree of the map (0 /a=1 , • • • , 4>' a=N2 ) —> (V\ a=1 , • • • , Vi =n2 ) seems some- 
what complicated. In order to get a more explicit form of Yn, let us move to the expression 
of (J)' by eigenvalues and SU (N) angles: 



<f>' = U 



\ 



A 



U\ UeSU(N). 



(B.5) 



A? 



/ 



Then, 



p-l N 



W) ab = 9 ? EE AfAj^" 1 *^ with ee y/NirffU)*. (B.6) 



=0 i,j=l 



Note that, from the completeness of the basis {£"}: X]^Li(^ a )«i(^ a )^ = j^^u^jk, ^tj satisfies 



N 2 



(B.7) 



a=l 
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Using this, one can see that each of = 1, • • • ,N) is an eigenvector of ^(0') 

whose corresponding eigenvalue is g P Y^Zo AfAip -1 . Hence, 

N / p-1 



t\ab 



i,j = l \ £=0 , 

^ n( W r)n(ff)' (b.8) 

The measure of 0' is expressed in terms of the eigenvalues and the angles as 

N 2 ,, la N 

= J]^ = ^(11^) A(A) 2 rff/, (B.9) 



, /2tt 

a=l v i=l 



where A(A) = rL>j(Aj — Aj) is the Vandermonde determinant, and dU is the SU(N) 
Haar measure normalized by J dU =1. Cn is a numerical constant depending only on 
iV. Then, the (//-integrals in (IB.3j) becomes 



A?" 2 

a=l V^^" 

. „ AT AT , 

= C w y dl7 y (J{d\) Hiw.-w^e^^^^, (B.10) 



j=l i=l i>j 

where Wi = g p X^ gives the Nicolai mapping. Thus we find 

N 2 

:n 



k!a \ ..2 



^ / du n 



1>J 

A' 



x]j5 ((tfn B U) u ) . (B.ll) 



i=l 



is the mapping degree defined by (1A.50I) . Plugging this into (1B.3|) . we end up with 

n(^-5(4;)n^(^)«) 

x E/ du ( un 'B Uj T (J d ^ J d $' a ft a ft a J ■ ( B - 12 ) 
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We changed the integration variable Q B as Q' B = U^VtgU under which the measure is 
invariant: dVt B = dQ' B . 

In (IB.12I) . the [/-integrals are clearly finite. The integration region for Q' B is the 
(iV 2 — l)-sphere S^ 2 " 1 defined by Yldj=i l(^s)ij| 2 = jj> an d the support of the delta- 
function n»=i ^ 1S a region of the S^ 2-1 determined by (Q B )u = ■ ■ ■ = (Q' b )nn = 
0, i.e. S^-^ -1 . The ^-integrals are also finite, because the integrand is a polynomial 
of (fijg)y multiplied by the delta-functions. Thus, the prefactors of the vanishing Grass- 
mann integrals in (1B.12j) are finite, meaning that Y/v definitely vanishes due to the trivial 
Grassmann integrals. 
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